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The focusing action of horizontal elongated irregularities of electron density in the 
ionosphere is considered. Expressions are developed for the general case where the azimuths 
of the radio star and of the short dimension of the irregularity may each be different from the 
azimuth of the interferometer baseline. These expressions appear to explain the main 
features of broadband radiostar scintillations, namely (a) their bandwidth, (b) their position- 
shift patterns, and (c) the nature of their association \\ it h spiead-F, To reconcile the theory 
with the occurrence data on scintillations which exhibit position shifts, it is necessary to 
postulate an irregularity movement towards an azimuth of 120°, i.e., a movement along 
the isoclines. 



1. Introduction 

Part I of this series [Singleton, 1964] discussed 
observations of Cassiopeia A made with the Boulder 
swept-frequency interferometer which operates in 
the frequency range 7.6 to 41 Mc/s. The scintilla- 
tions observed were found to be broadband, to 
involve occasionally position shifts and dispersion 
and to be associated closely with spread-F. The 
aim of this part is to investigate to what extent 
refractive properties of ionospheric irregularities 
can be invoked to explain these observations. 

2. Model Employed to Interpret 
Position Shifts 

The broadband nature of the scintillations first 
observed by Wild and Roberts [1956] in the fre- 
quency range 40 to 70 Mc/s led them to postulate 
that the ionospheric irregularities responsible be- 
haved as lenses. They assumed that the irregulari- 
ties were physically larger than the first Fresnel 
zone and that as a consequence the deviations pro- 
duced in an incoming plane wave front were due 
mainly to refraction rather than diffraction. In 
other words, t^ie ionospheric lenses were considered 
to fqcus the energy of the radio star so as to produce 
simultaneous increases of intensity over a large 
range of wavelengths. Warwick [1964], working 
in the frequency range 7.6 to 41 Mc/s, has recently 
shown that such a focusing mechanism is consistent, 
not only with the broadband nature of the scintilla- 
tions, but also with the configuration of position 
shifts which are occasionally associated with scin- 
tillations in this frequency range. 

i On sabbatical leave from the Physics Department, University of Queensland, 
Australia. 



Warwick's analysis involved the rather special 
case in which the line-of-sight to the source, the 
interferometer baseline and the short dimension of 
the focusing irregularity (assumed elongated and 
horizontal) are taken as coplanar. Here an analysis 
is outlined which is similar to Warwick's but which 
deals with the general case where the azimuth of the 
source and that of the short dimension of the ir- 
regularity may each be different from the azimuth 
of the interferometer baseline. In the first instance 
a flat earth approximation will be discussed for 
elongated horizontal irregularities. The effect on 
the results of changing to a spherical earth and 
ionosphere and allowing the irregularities to be 
other than horizontal will then be considered. 



3. Focusing Action of Horizontal Elongated 
Irregularities Above a Flat Earth 

By assuming that the horizontal elongated irreg- 
ularities are moving over a flat earth, we simplify 
the problem to the extent that all directions retain 
their identity when attention is transferred from 
the level of the irregularity to that of the observer. 
The three basic directions in the problem are the 
direction of the baseline reckoned from its southerly 
end to its northerly end, the direction of the radio 
source and the direction of movement of the focusing 
irregularity. At this stage it is assumed that the 
velocity of the irregularity is horizontal and normal 
to its larger dimension. These three directions are 
characterized by azimuths (measured from the north 
through east) a (baseline), f 5 (source) and v (velocity). 
For Casseopeia A as seen from Boulder, i- s lies 
between 317° and 43°, a is constant (41°), and v may 
have any value between 0° and 360°. 
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3. 1 . Geometry of the Refracted Ray 

In figure 1 let SO be a ray propagating in a medium 
of refractive index ju+5/x which encounters at a 
boundary separating this medium from one of the 
refractive index fi. It will be refracted at to 
propagate in the second medium along OP. Let 
OS' be the continuation of PO, then the angle 
through which the ray is refracted is f$ e . 6 S and % s 
represent the zenith angle and azimuth of the 
incident ray while 6 e and £ e represent the zenith 
angle and azimuth of the refracted ray. A con- 
sideration of the spherical triangle ZSS' yields the 
relationship 

cos (&>—&) = (cos ft,— cos 6 S cos e )/sin 6 S sin 6 e . (1) 

An adaptation of SnelFs law [Warwick, 1964] 
gives the approximate formula 



sin ft = (5/x//i) tan p 



(2) 



where p is the angle of incidence of the ray SO. 
If we neglect the earth's magnetic field, the refractive 
index is given by 

where /o is the plasma frequency and/ the frequency 
of the propagating wave. Thus (2) may be rewritten 
as 



sin^=-{/o5/o/(/ 2 -/S)}tanp. 



(3) 



Here (/o+<5/ ) is the plasma frequency in the medium 
in which the incident ray propagates and j the 
plasma frequency in the medium in which the 




Figure 1. Geometry of an incident and refracted ray SOP. 



refracted ray propagates. For the values of / 0r 
8f , and / commonly encountered viz, 5 Mc/s, 0.3 
Mc/s, and 30 Mc/s, Mfo/if-fo) is of the order of 
0.002. Thus for values of p up to 88.5°, ft is less 
than 2.5° and cos ft equals one with an error of less 
than one part in a thousand. This error increases 
to one part in a hundred for values of p between 
88.5° and 89°. Consequently, except for angles of 
incidence greater than about 89°, cos ft may be 
equated to one with little loss of accuracy. There- 
fore (1) may be written 

cos (£,— £ s ) = (l— cos 9 S cos 0<,)/sin 6 S sin B e . 

For the case where the radius of curvature is 
horizontal it can be shown that 

cos 6 e =ix cos d s /(n-\-8jj). 

This case is an accurate representation of the situa- 
tion for rays passing through the edge of an iono- 
spheric lens. Since 5/x//x is of the order of —0.002, 
/x/(/x+5/x) is of the order of 1.002 and 6 e is virtually 
equivalent to 6 S . Thus cos (£ e — f s ) is very close to 
unity. A detailed examination of (1 — cos d s cos 6 e )/ 
sin 6 S sin 6 e for the range of values of 6 S involved in 
the observations (50° to 81°) shows that cos (£ e — £ s ) 
is equivalent to unity with an accuracy of better 
than 0.1 percent. Thus little error occurs if changes 
in azimuth on reflection are neglected except when p 
is nearly 90°, that is when the irregularities are 
viewed along their length. 

3.2. Geometry of a Received Ray 

In figure 2 let OSQ be a vertical plane containing 
the received ray SO which has an azimuth £ and a 
zenith angle 6, OB is the direction of the interfer- 
ometer baseline which has an azimuth <r. Then 
^=£— <r. Let YO be perpendicular to OB in the 
plane SOB. Then, <p is the complement of the 
angle measured by the interferometer, i.e., y. 




Figure 2. Geometry of a received ray SO. 
OB is the baseline direction. 
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In the spherical triangle BSZ it can be shown that 
sin <p= cos 7= sin 6 cos \p. 

When all such received rays have a constant azimuth 
it follows that small changes in 6, involved in moving 
from one ray to another, are related to the corre- 
sponding small changes in <p by 



5<p=86 cos 6 cos i^/sin y. 



(4) 



3.3. Edge Rays 

In figure 3, DKL is the line image produced by the 
ionospheric lens whose long dimension is in the direc- 
tion of CH and JI and whose short dimension is in 
the direction JC. The interferometer baseline is 
represented by MN (azimuth a) and UFMK and 
WHNL are focused rays which pass through the 
interferometer antennas M and N. Since the 
direction of arrival varies monotonically across the 
baseline at a nearly uniform rate, the " average ray" 




Figure 3. Geometry of edge rays. 
The ionospheric lens JIIIC produces the line image DKL. 
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to whose direction the interferometer responds is GO 
where O and G are the midpoints of FH and MN, 
respectively. 

In section 3.1 it was shown that, except along the 
long dimension of an irregularity, the azimuths of 
rays from near the edge of the lens are unchanged 
from that of the source. Consequently FM, GO, 
and HN have the same azimuth as a ray which 
passes undeviated through the central portion of the 
lens (AX). However, the rays FM, GO, and HN do 
not necessarily meet at a point on the line image DL, 
since the azimuthal direction of the source and the 
baseline direction are, in general, skew. 

Rays BP, EQ, and CR are parallel to FM, GO, 
and HN respectively, in the vertical cut of the lens 
IACD which has the same azimuth as the source 
(J,). Rays BP, EQ, and CR are obtained from FM, 
GO, and HN by projecting along the direction of the 
long dimension of the lens. Then BP, EQ, and CR 
necessarily intersect on the line image in a point D 
and 6 S , O , and 6 e are thus the zenith angles of the 
source (i.e., the undeviated ray AD, EQ (or GO) 
and CR (or HN) respectively. Line PQR is the 
projection of the baseline on the azimuthal plane 
IACD, while XYZ is parallel to the short dimension 
of the irregularity JC (azimuth v) and is therefore 
normal to RN. 

Let SN be perpendicular to MN in the plane 
MNH and TO perpendicular to MO in the plane 
MOG. Thus <p e and <p are the complements of the 
angles which the interferometer would measure for 
the directions HN and GO. Let <p s be the comple- 
ment of the angle the interferometer would measure 
if the source were viewed directly in the absence of 
the lens. 

Considering the projection of MN on XR we have 

PR = MN cos (<r—p)/cos (v— f,) = A cos (^— k)/cos k 



where $=<?— £„ k=v~ £ s and A is the length of 
baseline. Detailed examination of this relationship 
for all values of % s and v shows that it is independent 
of the magnitude of v but changes sign when £ s is 
larger than 180°. Therefore, in general 



PR=±A COS {yp— it) /COS K 



(5) 



where the positive sign applies for £ 5 <180° an d the 
negative sign applies when £ s ^>180°. 

Considering the edge ray HN and the average 
ray GO, it can be seen that 

(<p e —(Ps) — (<fo— <Ps) = <Pe~ <P0- 

Since <p e is very little different from <p s or <p , use can 
be made of (4) to produce the following: 

(8 e —6 s ) cos 6 S cos \f//sm y tS —(fpo—<ps) 

= (6 e —do) cos d s cos ^/sin y S) (6) 

6 e —6 s , <PQ—(p s and 6 e — 9 are functions of frequency 
d e — 9 S will be written 80 e (f) and <pq—<p s as 8<p (J) or 
— 8y (f). 8y Q (f) is the change in the apparent 
position of the source, as measured by the interfer- 
ometer, as the edge of the irregularity just comes into 
view; 6 e —d can be determined as a function of 
frequency as follows. 

In figures 4a and 4b we have drawn the vertical 
plane through the source for observing frequencies 
( /) larger than the focus frequency (f m ) and smaller 
than the focus frequency. The focus frequency is 
defined [Warwick, 1964 J as the frequency for which 
the focal length of the ionospheric lens [F(f)] equals 
its range (d). It can be seen from figures 4a and 4b 
that 

9-d =±Pn/2{F(f)-d} 





(0) 



(b) 



f < K 



Figure 4. Illustrating the relationship between 6 e and O and the focal length F(f) and 
range d of an irregularity for (a) f>f m and (b) f<C/ m . 



1098 



where the positive sign applies to the case/>/ w , and 
the negative sign to the case f<Cfm> In what follows 
only the case />/ m will be considered. With (5), 
(6) can now be written 



86 e {f) cos 9, cos ^/sin y 8 +ho(f) 



■-± 



COS (\f/—K) cos 6 S cos ^ 



2{^C/)-rf} cos k 



sin 7, 



(7) 



Since the refraction of edge rays results in changes 
in zenith angle rather than azimuth (sec. 3.1), and 
even these changes are small, (3) may be written 



or 



W.(/) = -{Mfo/(f 2 -f 2 o)} tan P 
8d e (f)ocll(f-f ). 



It follows that 



F{j)oz(f-fl) 



(8) 



where i*\/) is the focal length of the lens. However, 
for most cases, /is an order of magnitude larger than 
/o and thus 

F(f)=df/fi (9) 

where d is the focal length and the range of the ir- 
regularity at the focus frequency f m . 

The angle measured by the interferometer y is 
given by 

cos y=(nc/J+p)/A 

where n is the fringe order, c the velocity of light, 
/the observing frequency, p the difference in length 
of the two antenna feeders and A the length of the 
baseline. It follows that 



8y=nc8f/f 2 A sin y 



(10) 



where 57 is the shift in the measured angle corre- 
sponding to an observed shift of 8/ in the frequency 
at which a given fringe order is excited. Since 
7i=// A/' where A/ is the fringe spacing, (10) can be 
written 

8y=(c/jAsmy)(8flAf). (11) 

Using (9) and (11), it is possible to rewrite (7) as 
follows 

ffn 



1 -7-,= a 



+ b 



where 



and 



a=± 



A/ " (/ 2 -/ 2 J ' 

A 2 cos (\//—k) cos 6 s cos \p 
2cd cos k 



-(A/c)f86 e (f) cos 0, cos ^. 



(12) 



(13) 



The positive si^n in the expression for a applies if 
£ s <180° and the negative if k>180°. Since 8d e (f) is 
proportional to/ 2 it follows that J" 2 dd e (f) is constant 
and so all the quantities involved in the expressions 

for a and b are constants for a particular scintillation. 
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Figure 5. Scatter diagram .showing the connection between 
f«f/Af and P f 2 ,n/(f 2 -'f 2 M,) for the scintillation at 0710 UT, 

1>, December, 1960. 



Thus fifl A/ should be a linear function oifimlU 2 — 
/„})• r Phis has been verified for a number of scintil- 
lations. Figure 5 shows the plot of fdf/Af versus 
ffm!(f—fm) for the scintillation at 0710 hr UT on 
14 December 1960 which was discussed in detail by 
Warwick [1964], The value of f m used (28 Mc/s) 
was obtained by an interpolation between the fringe 
slopes on either side of the fringe-slope changeover, 
relative to the star's diurnal motion. It will be 
noted that tin 1 expression for a, (12) contains the 
range of the irregularity d. Unfortunately evalua- 
tion of a does not lead immediately to d since the 
direction of movement of the irregularity and hence 
k are unknown. 

The following consideration of figure 6 will yield 
expressions for the duration of the scintillation and 
the width and velocity of the focusing irregularity. 
The duration t(j) of a scintillation at frequency /is 
given by 

t (/)=/? cos */« 

where * is the irregularity, and hence the image, 
velocity, D is giveYi by 

z?=± h F(f)-d }//''(./)= ± 1 (f-f m yfi 

where the positive sign applies to the case />/» 
(fig. 6) and the negative sign to the case f<CJm- 
Thus ' 

t(f) = ±l cos K<J*-f m )l~P. 

As/ increases (in the/>/ w case), t(f) reaches a limit- 
ins: value t, and in terms of this 



t(j)=±t i (f-f m )/f. 



(14) 



Thus as frequency increases t(f) decreases from ex- 
tremely large values for small values of / to zero 
at f m and increases again to t t for large values of/. 
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734-668—64- 




Figure 6. This figure illustrates the derivation of expressions 
for speed of movement (-») and irregularity width (I). 



This offers an explanation for the narrowing of the 
scintillation near/ m as reported in part I. 

From figure 6 it is evident that the breadth of the 
focusing irregularity (i) is given by 

L=2 OR=2 OP cos k/cos 6 8 =2F(J)80 e (J) cos /c/cos d s 



i.e., 



L=2d86 e (f m ) cos /c/cos 6 S 



(15) 



Making use of (13) this becomes 

L=—2bcd cos \i\p m A cos 2 6 S cos \p. (16) 

All the quantities but k on the righthand side of 
this equation are obtainable for any scintillation and 
thus L cannot be calculated without a knowledge of 
the direction of movement of the irregularity. 

The magnitude of the irregularity velocity is ob- 
tained directly from a knowledge of L and t t since 

v=L/ti=—2bcd cos K/tifmA cos 2 d s cos \p. (17) 

Again a lack of knowledge of k prevents evaluation 

of v. 

3.4. Apparent Angular Drift of the Source 

As a focusing irregularity moves across the line 
of sight of the radio star it will change the apparent 
position of the star in a definite pattern, as explained 
by Warwick [1964]. In this section an expression 
will be derived for the rate of change of the observed 
angular position of the source. At the beginning of 
a scintillation the observed angular position of the 
source changes abruptly from y s to y (f) (for the 



edge ray) that is by an amount 5y (f). There is 
then a gradual recovery to y s and then a gradual 
overshoot to y s — 8y (f). Thus the total change in 
the apparent angular position of the source is 
25y (f) and this occurs in a time t(f). The rate of 
change of the observed angular position of the 
source W (f) is consequently given by 

Wo(j) =2dy ( f)lt(f). 

From (7), (9), (14), and (15), W (f) can be written 
finally as 



Wo(f)- 



n Jm 



(J 2 -f m ) 
r 



{ 



} 



I 



cos 6 S cos (£ s — a) 
d cos y s 

A f 
-COS 6 S ± T W ri cos (v—c 

Lf-Pm _ \ (18) 

cos (£ s — v) 



i}<„ 



where the positive sign before the second term applies 
if £ s is less than 180° and the negative if % s is larger 
than 180°. This expression applies for j^>j m but 
can be converted into the analogous expression for the 
case f<Cfm by reversing the signs before the second 
term. 

At frequencies higher than the focus frequency, 
the CB H and CA H type scintillations involve values of 
W (f) which are of opposite sign. It is instructive, 
therefore, to examine the variations of W (f) with 
source azimuth £ 5 and wind azimuth v. In figure 7 
the full line and dashed curves represent W (f) 
(J 2 —jm)l fm plotted against source azimuth for 
several values of the wind azimuth. Here/ has been 
taken as 1.2/ m , L given the value 5 km, and only the 
range of azimuths and elevations appropriate to the 
observations of part I have been considered. The 
curves have an infinity when £ s =v— 90 c . Under this 
condition the irregularity is being viewed along its 
long dimension and its focusing action will be such 
as to produce changes in azimuth rather than ele- 
vation. However, the approximation on which the 
curves of figure 7 are based assumes that the focusing 
action involves changes in elevation only (sec. 3.1). 
Thus the infinities in the curves at % s =v — 90° are 
unlikely to be meaningful. 

To interpret correctly the curves of figure 7 it is 
essential to know over which parts of the curves the 
approximation on which they are based is accurate. 
An accuracy criterion for these curves can be estab- 
lished in the case where the normal to the refracting 
surface, where an edge ray strikes it, is horizontal. 
In this case it can be shown readily that 



cos p=sin d s cos (*>—&). 



(19) 



It was pointed out in section 3.1 that it is reasonable 
to assume that there is negligible change in azimuth 
on refraction when cos j3 e is approximately equal to 
one, that is when p is less than 88.5°. Using this 
limitation on p it is possible to employ (19) to find 
the range of values of £ s , for each v, for which the 
approximation used in the construction of figure 7 is 
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Figure 7. Variation of the rate of change of the observed 
angular position of a refracted ray with source azimuth (£ H ) 
for several values of wind azimuth (v). 

accurate (full lines). As anticipated, the approxi- 
mation, based on no change in azimuth on refraction, 
breaks down (dashed lines) in the vicinity of the 
infinities at which the sign of W (j)(f 2 -f 2 m )/ f 2 m 
changes. Thus, in those cases where W (f) (j 2 —jm)l 
t<j 2 m changes sign it is only possible to determine 
part of the curve which represents the variation of 
this quantity with £,. A complete theory would, of 
course, allow all the curve to be determined. Lack- 
ing such a theory, we assumed that a connection 
between the two branches of the curve, consistent 
with the curvature trends of the branches, approxi- 
mates to the true variation. These connections are 
drawn as dotted lines on figure 7. 
It can be shown that 
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Assuming that the smallest detectable value of 8f is 
±0.1 Mc/s, we can determine how the smallest 
detectable value of ^o(/)(/ 2 -/m)/^ varies with £ s . 
The dash-dotted lines of figure 7 represent this 
variation for an/ m of 30 Mc/s and an / of 36 Mc/s. 
In a CB H [Singleton, 1964] scintillation the fringe 
drift is from high to low frequency as time progresses. 
Now 

cos yoGO =nc/fA+p/A. 



I f n is kept constant and nonzero and / allowed to 
decrease, obviously y (f) decreases. But W (j) is 
equivalent to <fyo(j)/dt so that CB H scintillations are 
associated with negative values of W (f). A similar 
argument shows that CA H scintillations are associated 
with positive values of M '„( /'). From diagrams such 
as figure 7 it is possible, therefore, to determine the 
ranges of values of source azimuth (£,.) over which 
detectable CB H and CA H scintillations will be pro- 
duced by a given wind azimuth (V). Scintillations 
occurring at those values of £, and v for which 
W (f) is below the limit of detection will correspond 
to type N as defined in part I. A summary of the 
ranges of values of £, over which CB H , CA H , and N 
type scintillations would be expected for a full 
range of values of v is presented in the lower part of 
figure 8. The upper two histograms in figure 8 
show how the occurrence of CB H and CA H scintil- 
lations vary with source azimuth. They show that 
there is a changeover from CB H to CA H at a source 
azimuth of about 8°. Such a change is predicted 
most accurately by a wind azimuth of 120°. 

In order to gain more information about the direc- 
tion of the wind, an attempt has been made to re- 
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Figure 8. The theoretical predictions of position shift proper- 
ties of scintillations {lower figure) are here compared with 
the experimentally determined histograms of CBh and CA H 
occurrence {upper histograms) . 
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Figure 9. Predicted histograms (broken lines) of (a) CB H , 
(b) CA H and (c) N type scintillations compared with the 
experimental histograms {full lines). 
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Figure 10. A polar plot of the relative occurrence of wind 
azimuths. 



construct the experimental histograms of CB H and 
CA H occurrence versus sidereal time. This was done 
by allotting weights to the predictions of figure 8 for 
wind azimuths of 20° to 160°. The weights^ (num- 
bers in brackets on fig. 8) were chosen by trial and 
error so as to get reasonable agreement between pre- 
diction and experiment for the CB H and CA H histo- 
grams. Figures 9a and 9b compare the experimental 
(full line) and predicted (broken line) histograms 
for CB H and CA H scintillations. 



Using the weights employed in the construction of 
figure 8, we can gain some idea of the variation of the 
relative occurrence with wind azimuth (fig. 10). 
Figure 10 also shows the azimuth of the isoclines in 
the observing region. It would appear that irregu- 
larities producing CB H and CA H type scintillations 
are usually elongated in the magnetic north-south 
direction, the component of motion normal to their 
length lying towards the southeast along the isoclines. 

Keverting to figures 8 and 9, we expect that the 
adopted weighting procedure would also reproduce 
the observed sidereal time variation of occurrence of 
the N condition. Figure 9c shows that this is not 
the case. However, it is interesting to note that the 
synthetic histogram (broken line), though smaller 
and narrower than the observed (full line) histogram 
reproduces, to some degree, the trends of the ob- 
served histogram. Note especially the slower rise 
of occurrence than fall off of occurrence with in- 
creasing sidereal time. From (20) it will be seen that 
a decrease of f m by a factor of three increases the 
minimum detectable value of W {j) (f—j 2 m)h' c fm by 
an order of magnitude. Such a change increases the 
range of £ s over which N type scintillations may be 
observed at the expense of CB H and CA H scintil- 
lations. Thus the synthetic histogram of figure 9c 
may be widened and increased in height to fit the 
observed histogram more closely if it is assumed that 
values of j m lower than 30 Mc/s predominate. The 
relative heights of the synthetic and observed histo- 
grams suggest that the number of occasions when 
/ m <30 Mc/s to the number when / m >30 Mc/s 
should be of the order of 2 to 1. 

4. Effect of a Spherical Earth and Iono- 
sphere and Nonhorizontal Irregularities 

The above discussions have been based on the 
approximation that the earth and ionosphere are 
flat. Little change in the form of the analysis is 
necessary when a change is made to a spherical earth 
and ionospheric situation. The main difficulty that 
arises is that the azimuth of a vector at a point other 
than on the observer's zenith cannot be transferred 
without change to the observing point. This is due 
to the convergence of meridians. Computations 
have been made which allow the assessment of the 
error involved in assuming that the apparent azimuth 
(as seen from the ground) of a vector in the ionosphere 
equals its true azimuth in the ionosphere. These 
computations are summarized in figure 1 1 where the 
error for several true azimuths is plotted against 
sidereal time. An irregularity height of 300 km has 
been assumed. It will be seen that the error never 
exceeds 3.2°. Considering the crude method used 
to determine the wind direction in section 3.4, an 
error of this magnitude is probably of little 
significance. 

So far the irregularities have been assumed to be 
elongated and horizontal and are thus similar to those 
thought to be responsible for some types of midlati- 
tude spread-^ [Bowman, 1 959]. There is a considera- 
ble bod of evidence however, which points to the 
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Figure 11. Error involved in assuming that the apparent azimuth of a vector at a dis- 
tant point in the ionosphere is equal to its true azimuth at that point, plotted against 
sidereal time. 



existence of large field-alined 
latitudes [Peterson et al., 1 
Knecht and Russell, 1962]. 
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the dip is of the order of 75° is 
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irregularities at these 
955; Dearden, 1962; 
Such irregularities, if 
inns, would appear to 
Since a region where 
being viewed at zenith 
to 80°, the geometrical 
r to that already con- 
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lodels with the experi- 
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5. Intensity of Focused Scintillations 

In the vicinity of the focus the beam of radiation 
focused by an ionospheric irregularity will have a 
cross section smaller than that of the irregularity. 
As a result the signal strength in the focused beam, 
near the focus, will be higher than the signal strength 
at the same point in the absence of the focusing 
irregularity. The gain in signal strength as the 
result of focusing can be calculated as follows. 

Consider figures 4a and 4b. The radiation which 
passes through an area of breadth AB reaches the 
ground over an area of the same length as the irregu- 
larity but with a breadth CD. The gain (G) in 
signal strength is obviously: 



G= 



CD 



1 



F(J) 

FU)~d 



(21) 



\vhere T the positive sign applies to the case where /> 
f m and the negative to the case where /</ m . Use 
(8) to evaluate F{J) ; then (21) becomes 

G=(i-/,v/ 2 )/ii-yy/ 2 i. 

In figure 12 the gain is plotted as a function of 
fjj for several values of f /f. The broken line parts 
of the curves represent the situation where f m is not 
greater than the critical frequency of the F-layer 
(fo). This is physically unlikely. Note the limited 
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Figure 12. Gain due to the focusing action of an ionospheric 
lens plotted against f m /f for several values of f /f . 
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range, about f m , of observing frequency (J) for which 
the gain is greater than one. Bandwidths up to 
7 : 1 are possible when the F-layer critical frequency 
/o is less than 10 percent of the observing frequency/. 
This upper limit to the bandwidth drops to 1.6 to 1 
for/o equal to 75 percent of/. The observed band- 
widths (part I) of 2:1 and occasionally 3:1 are 
consistent with this result. 

In figure 13,f m /fis plotted against f /j for a gain of 
one. Gains greater than one occur only for values of 
ijj and / // which correspond to points within this 
curve. On the left and bottom of the figure several 
scales of observing frequency / have been added for 
particular values of Jrn and f . The effective band of 
the swept-frequency interferometer (after the MUFE 
has been taken into account) is shown on these scales 
by the thick lines. It is evident that some or all of 
a scintillation whose focus frequency j m fulfills the 
condition 2 Mc/s>/ w >56 Mc/s may be observed 
with the interferometer. 



Figure 13 offers an explanation for those scintilla- 
tions which exist only in the low-frequency octaves. 
Consider a focus frequency of 5 Mc/s and an i^-layer 
critical frequency of 2 Mc/s. As the observing 
frequency (J) is allowed to vary from 10 to 40 Mc/s, 
curve (a) is traced out. It can be seen that under 
these conditions the gain is greater than one only for 
j less than 25 Mc/s. For the same value of j m but 
a higher value of f , say 4 Mc/s, curve (b) shows that, 
only for observing frequencies less than 20 Mc/s, 
the gain will be greater than one. Thus low- 
frequency scintillations apparently occur when j m 
and/o are small. Low values of / will be associated 
with the low values of MUFE necessary to permit 
observation of low-frequency scintillations. How- 
ever, as found experimentally in part I, low values of 
MUFE, while a necessary condition for the observa- 
tion of low-frequency scintillations, is not a sufficient 
condition. The focusing frequency must also be 
small. 
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Figure 13. Here i m ji is plotted against f /f for a focusing-action gain of one (full-line curve). 

Several /scales for particular values of f m and/ are also included. The broken-line curves show the variation of/ m // with /<,// for (a) f m = 5 Mc/s and f =2 Mc/s, 

and (b) f m = 5 Me/s and / G =4 Mc/s as / varies from 10 to 40 Mc/s. 
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6. Association Between Focusing Frequency 
and Spread-F Parameters 

The focusing properties of an ionospheric irregu- 
larity depend on its incremental electron density, its 
size and the nature of the curvature of its surfaces. 
By assuming that the irregularities which produce 
scintillations are the same as (hose which simultane- 
ously produce spread-^, it is possible to gain infor- 
mation about the incremental electron density of the 
scintillation-producing irregularities. 

Conventional ray optics will now be used to 
develop connections between the focusing frequency 
of an irregularity and its incremental electron 
density. Three curvature models will be considered. 
In the first model the irregularities will be assumed 
to be semicylindrical deviations from isoionic con- 
tours that are otherwise flat. Cylindrical columns 
will be considered as the second model. In the 
third model the isoionic contours will be assumed 
to have a sinusoidal wavelike nature. 

6.1. Semicylindrical Irregularities 

Figure 14a depicts an irregularity which takes the 
form of a semicylindrical deviation of an otherwise 
flat isoionic contour, it will be assumed that the 
refractive index immediately beneath the contour is 
fi, while that immediately above is /z + <5,u. The 
irregularity, which lias a radius of curvature r, i> at 
a distance d from the observing point. It is well 
known [Longhurst, 1957] that the focal length (F) 
of such a system is given by 



F r 



(22) 



Neglecting the effect of the earth's magnetic 
field, the refractive index of a plasma for an observing 
frequency/ is given by 



and 



M 2 =l-/o7/ 2 

(/x+5 M ) 2 =l-C/o+5/o) 2 // 2 



(23a) 
(23b) 



where /o is the plasma frequency below the boundary 
in figure 4a and /o+5/o is the plasma frequency 
above this boundary. Using these two expressions 
in connection with (22) and assuming that the 
radius of curvature is at least an order of magnitude 
smaller than the focal length, we find that 



(rlF)(f-fl)=28f f +Sf o . 



(24) 



The variation of plasma frequency at the peak of the 
i^-layer associated with most spread-i^ 7 configurations 
is an order of magnitude less than the minimum 
plasma frequency here. Thus (24) can be written 

8jol.fo=(r/2F){(My-l}- 

When the observing frequency/ equals the focus 
frequency j m the focal length F equals the range d 
of the focusing irregularity. Thus 



«/o//o=0'/2r/){(y m //o) 2 -l} 



(25) 



which is the required connection between the focus- 
ing frequency and the parameters of the associated 
spread-^ configuration. 

If, as in figure 14a, the irregularity is convex as 
viewed from below, d and r are required to have 






(a) 



(b) (c) 

Figuke 14. The three irregularity models whose focusing properties have been investigated. 
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opposite sign (according to the sign convention used) 
to enable the lens to be convergent and thus 5/ is 
negative. If the irregularity is concave as viewed 
from below <5/ has to be positive for the lens to be 
convergent. 

6.2. Cylindrical Irregularities 

The focusing action of a cylindrical irregularity in 
a plane normal to its axis will now be considered. 
Assuming that the irregularity has a circular cross 
section of radius r and that the refractive index of 
the plasma within the irregularity is m + ^m and that 
outside m (fig. 14b), simple ray optics [Longhurst, 
1957] can be used to show that the focal length F 
is given by 

v 6/jl 

where F and r have to be of opposite sign for the lens 
to be convergent. Since F is expected to be at least 
an order of magnitude larger than r, this may be 
rewritten 



2F/r=n/8», 



(26) 



which is of similar form to (22). Employing the 
methods outlined in section 6.1, but using (26) 
rather than (22), we obtain the following connection 
between the focusing frequency, range of the ir- 
regularity and the parameters of the associated 
spread-F. 



f=-^{(^)-i: 



(27) 



Here d and r have to be of opposite sign for the lens 
to be convergent, and hence 8f is negative. 

6.3. Sinusoidal Irregularity 

Warwick [1964] considered the focusing action of 
irregularities resulting from sinusoidal wavelike 
isoelectronic contours. Such irregularities are de- 
picted in figure 14c. He demonstrated that rays 
passing through the lens at points distant \y\ from its 
axis are brought to a focus at a distance / given by 

f=(\/2TL)y/(sm2wy/\) 

where L=8eA /2e and A , A, and the dielectric con- 
stants e and e + de are as shown in figure 14c. 
Because of the caustic effect associated with such 
lenses, Warwick asserts that the rays passing through 
the edge of the circle of least confusion arise from 
the points y=±\/io, and the effective focal length is 
thus 



Since 



5.878ttA <5e 
■A-fl/f, it follows that 



8ele=-(28j /f Q )/{(f/f o y- 



and thus 



_x* (/ 2 /io)-i 

5.878ttA 25/ //o 



At the focus frequency (f n ), F equals the range of 
the irregularity (d). Consequently, the final con- 
nection between the focus frequency, the range and 
dimensions of the irregularity, and the parameters of 
the associated spread--? 1 configuration is 



5/0. 



/o 36.9.4 d 



{©-}• 



(28) 



6.4. Experimental Verification 



In sections 6.1, 6.2, and 6.3 it has been demon- 
strated that the relationship between the focus 
frequency (1 m ), the range of the irregularity (d) and 
the parameters of the associated spread-.? 1 con- 
figuration (/ and <5/o) is 

Mitt)'-} 

where 6' is a function of the size and shape of the 
irregularity. In part I, figure 15, it was shown that 
/o<5/o is proportional to j 2 m . This is consistent with 
(29) since j m is usually an order of magnitude larger 
than / . r lhe same body of data in part I, figure 15, 
is used in a scatter plot (fig. 15) of 8f /J versus 
{ (fm/jo) 2 — 1 }• A linear relationship obviously exists, 
the correlation coefficient of the scatter being 0.79 
with a 1 percent confidence level of 0.41. The slope 
of the line fitted to the scatter is 8.57X10" 4 . 

Irregularities at 400 km viewed between 1400 and 
1 700 hr sidereal time are at a range (d) of 930 km on 
the average. It follows that 6' equals 0.80 km on the 
average. Values of the breadth of the average 
focusing irregularity can be calculated in terms of s. 
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These calculations are summarized in table 1 . En the 
case of the sinusoidal irregularity two widths are 
quoted. The first of these is calculated on the basis 
that the wavelength of the isoelectronic contour is 
comparable with its amplitude as assumed by 
Warwick [1964]. The second width quoted is based 
on the assumption that the depth of the focusing part 
of the irregularity is half its width, a situation which 
is directly comparable with the case of the semi- 
cylindrical irregularity. 

The radius of the first Fresnel zone at a range of 
930 km for the wavelengths used in most of the 
observations is 2.5 to 3.5 km. Thus it appears that 
the irregularities must be either of the cylindrical or 
sinusoidal types in order to be larger than this zone 
and allow refraction effects to dominate possible 
diffraction effects. 

Table 1. Comparison of the irregularity width for the three 
models proposed 



Semi cylindrical 
irregularity 


Cylindrical 

irregularity 


Sinusoidal 

incuularity 


,S=r/2=0.80kni 

r=l.Gkm 

Width of irreg. X=2r=3.2 
km 


,S= r /4 = 0.8() km 

r = ::5.2kin 

Width of Lrreg. L=2r= 
6.4 km 


S=X2/3(UKlo = 0.8<) km 

XV< 1 = 29.5 

Width of lrreg. L=X/2= 
14.7 km (.1„=X) 
= X/2=3.7 km (,1„ = 
A/4) 



7. Dispersion 

As Wild and Roberts [1956] first suggested, the 
dispersion effects noted in part I are apparent ly 
due to horizontal gradients in electron density which 
act as huge prisms. Since scintillations with dis- 
persion exhibit much the same position-shift prop- 
erties as other broadband scintillations, focusing 
structures of the types discussed above as well as 
horizontal gradients are required for their explana- 
tion. 

No complete analysis of dispersion by horizontal 
gradients of electron density will be given here. 
Attention is drawn, however, to a possible explana- 
tion of the sidereal time A r ariation of the dispersion 
effects. It was noted in part I that scintillations 
are frequently associated with dispersion before 
0900 hr sidereal time but that the incidence is low 
from this time until 1600 hr when observations cease. 
Figure 16 shows the track of the hne of sight to 
Casseopeia A across the 200, 400, and 600 km 
levels. Each track has been divided up into 2-hr 
intervals of sidereal time. The figure also shows 
the 70°, 75°, and 80° isoclines (lines of constant 
magnetic dip). It is obvious from the figure that 
an isocline corresponding to a dip of 77° would 
cross the ^=600 km track at about 0900 and 1600 
hr sidereal time, whereas the 75° isocline crosses the 
h =400 km track at these times. The isocline 
corresponding to a dip in the vicinity of 75° may 
therefore form the northern boundary of the region 
with horizontal gradients of electron density suitable 
for the production of dispersion. 




105° 
West Longitude 



Figure 10. This shows the track of the line-of-sight from 
Boulder to Casseopeia A across the 200, 400, and 600 km levels 
(full lines). 
Each of the (racks have been divided up into 2-hr intervals of sidereal time. 

The broken linos represent the 70, 75, and 80° isoclines and the places marked are 

i he local ion of ionosondes. 



8. Conclusions 

The focusing action of large scale elongated 
ionospheric irregularities has been considered for the 
genera] case where the azimuth of the radio source 
and that of the short dimension of the ionospheric 
irregularity may both be different from that of the 
interferometer baseline. r l ne expressions developed 
appear to explain the main features of the broadband 
scintillations reported in part I, namely, (a) their 
broad bandwidth, (b) their position shift patterns 
and (c) their association with spread-/ 1 . To re- 
concile the theory with the statistical information 
concerning position shifts it was found necessary to 
postulate an irregularity movement towards an 
azimuth of 120°. This corresponds to movement 
along the isoclines with the long dimension of the 
irregularities in the magnetic north-south direction. 
The results do not permit differentiation between a 
horizontal irregularity model and a field-alined 
irregularity model. 

Because of a lack of knowledge of the direction of 
motion of individual scintillations it is not possible 
to find directly the range, width, or speed of the 
irregularities. 
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